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ABSTRACT: We continue to study the holographic QCD (hQCD) model, proposed in a 
previous paper, in an Einstein-Maxwell-Dilaton (EMD) system. In this paper we discuss 
some aspects of quark gluon plasma (QGP) in the hQCD model, such as drag force, jet 
quenching parameter and screening length. The results turn out to be consistent with 
those as expected in QCD qualitatively. By calculating free energy of the background 
black hole solution, we find that there exists a phase transition between small black hole 
and large black hole when chemical potential u is less than the critical one ue, and the 
phase transition is absent when chemical potential is beyond the critical one. 
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1 Introduction 


The analysis of various dynamical quantities from the experimental data obtained at Rel- 
ativistic Heavy Ion Collider (RHIC) in Brookhaven National Laboratory leads to strong 
evidence for a strong coupled quark gluon plasma (QGP), a deconfined phase of QCD at 
high temperature and high number density [1], [2], [3]. In the Au+Au collision with the 
maximum center-of-mass energy around 200 GeV, several phenomenological features, e.g, 
a very small value of shear viscosity, quenching of high energy partons with large transverse 
momentum, and elliptic flow etc., indicate that the dynamics of thermal medium produced 
after collision is dominated by non-perturbative effects [4]. Being formulated in Euclidean 
time, lattice QCD seems to be the best candidate to explain the phenomena in thermal 
equilibrium, while the perturbative QCD works only in the weak coupling region. But, 
both the perturbative QCD and lattice QCD are failed to compute some dynamical quan- 
tities like transport coefficients, drag force, and jet quenching parameter etc., in the strong 
coupling regime. Thanks to the feature of strong/weak coupling duality, the AdS/CFT 
correspondence [5] [6] provides a powerful tool to stress those issues. 

In the AdS/CFT correspondence, a well-known example of strong/weak duality is 
AdS;/CFT, and in this case, the four dimensional conformal field theory is M = 4 SU(N.) 
super Yang-Mills (SYM), while the bulk theory is the type IIB supergravity in AdS5. In 
the AdS;/CFT, correspondence, the ratio of shear viscosity over entropy density for the 
dual field theory is found to be small (ņ/s = 1/47) for large number of colors (Nc) and 


large ’t Hooft coupling (A = gy m° No) [7] [8] [9]. It turns out that this value is consistent 
with the one from RHIC data with experimental error [10]. In addition, it is found that the 
value 1/47 is universal for various deformations of the gauge theories with gravity duals. 
Another example of the universality is discussed in [11], where it is found that for some 
gauge theories (not necessarily conformal) with well defined gravity duals where if some 
conditions on the bulk stress energy tensor are satisfied, the electrical conductivity at finite 
chemical potential (u) and temperature (T), the thermal conductivity (xr) and the ratio 
of thermal conductivity to viscosity (Fu) are independent of any specific model. 

Certainly a realistic holographic model dual to the strong coupled QCD at finite tem- 
perature and finite density is a good starting point to study the RHIC physics. Unfor- 
tunately such a model is still absent. However, due to the universality mentioned above, 
one is expecting to understand various features of QGP at RHIC with some deformed 
AdSs5/CFT4, such as drag force, jet-quenching parameter, and screening-length etc. In- 
deed, over the past years a lot of works have been done along this line. 

By employing the so-called potential reconstruction approach [12-15] to build up holo- 
graphic QCD model [16] from bottom up point of view, a hQCD model in an Einstein- 
Maxwell-dilaton theory is proposed in [15]. In [14, 15], some properties of the hQCD 
model are studied, such as equation of state, Wilson line operators [12] and confine- 
ment/deconfinement phase transition and associated phase diagram [15], with good agree- 
ment with lattice QCD simulation. It shows that this model can capture some most impor- 
tant characteristics of realistic QCD. In this paper, we continue to study the hQCD model 
by investigating some aspects of QGP phase such as drag force, jet quenching parameter as 
well as screening length. The results show that these quantities are consistent with QGP 
properties. In this sense, we further confirm that this model may give us the hints on the 
efficient way to deform the pure AdS; geometry and to realize the holographic description 
of low energy QCD. 

When high-energy partons perform a dragged motion as they pass through the QGP 
medium, their energy loss can be encoded by drag force. In the AdS/CFT correspondence, 
an external heavy quark and its gluonic neighborhood are mapped into the endpoint of 
fundamental string attached to the AdS boundary and the string itself in the AdS bulk 
geometry, respectively. This external quark, with mass proportional to the length of the 
string, loses its energy as the string trailing back imposes a drag force on it. Within the 
framework of gauge/gravity duality, the drag force experienced by an external heavy quark 
moving with a constant velocity in M = 4 super Yang-Mills plasma at finite temperature 
is computed [17], [18], [19], [20], [21], [22], [23]. There are also further generalizations of 
drag force computation for the charged M = 4 SYM [24] and with the backreaction effect 
due to the static heavy quark cloud distribution [25]. 

Due to the medium, the suppression of heavy quark with high transverse momentum 
leads to energy loss which is so called jet quenching phenomenon [26], [27], [28]. The 
transport coefficient ¢ [29], characterizing such phenomena, is defined perturbatively as 
the ratio of square of the mean transverse momentum over the mean free path [30]. With 
the framework of eikonal approximation [31], [32], the parameter ĝ can be calculated from 
light-like Wilson loop in adjoint representation [33], [34], [35], [36], [37], [38], [39]. Ref. [40] 


and [41] calculate the jet quenching parameter in the presence of chemical potential. Using 
AdS/CFT, the quark-antiquark pair is mapped to the two endpoints of a fundamental 
string, both endpoints attach on the AdS boundary. The string configuration shows it hangs 
down to the bulk along radial direction and turns back to the boundary again. The light-like 
Wilson loop in the fundamental representation of boundary gauge theory and its thermal 
expectation value correspond to the trajectory of two endpoints of the string and exp(S) 
respectively, where S' is the string Nambu-Goto action. In addition, we further extend 
our analysis by studying the time-like Wilson loop and relate it with another important 
parameter, the screening length of a quark-antiquark pair. It is defined as the maximum 
length between a moving gg pair, beyond which they break off with no binding energy 
and thus become screened in the QGP medium. Screening length depends on the velocity 
and the orientation of the quark-antiquark pair with respect to the medium. Treating 
N = 4 SYM at finite temperature as a boundary theory, the binding energy between quark 
and antiquark pair moving in the QGP and the screening length are calculated in [42]. 
There are also other important generalizations to calculate these quantities by introducing 
a boost on the original static background [43], [44], [45]. In this paper we will calculate 
the screening length in the static frame of gq pair for our hQCD model. 

In [15] we calculated heavy quark potential between a quark-antiquark pair in our 
model, and found that there is a confinement/deconfinement phase transition, and that 
there is a critical point in the T — u phase diagram. In this paper we further confirm 
this phase transition in the small u region by computing free energy of the background 
black hole solution by using the method in [46]. In the AdS/CFT correspondence, the 
Hawking-Page phase transition [47] between AdS black hole and thermal gas in AdS space 
is identified with the confinement/deconfiement transition in gauge theory [48][49][50][51]. 
In our case, the phase transition happens between small black hole and large black hole, 
which will be clear shortly. In addition, we will argue that there is no phase transition 
between black hole solution and thermal gas solution in our model. 

The organization of the paper is as follows. In section 2 we briefly review the poten- 
tial reconstruction approach to the Einstein-Maxwell-Dilaton system by generalizing the 
discussion in [15] to the case with a coupling between dilaton field and Maxwell field. In 
section 3, we discuss the generic black hole solutions with asymptotical AdS boundary, and 
in particular present an analytic black hole solution. In addition, in this section we also 
briefly review the black hole solution for the hQCD model studied in [15]. In section 4, we 
calculate the drag force in this hQCD model. The jet quenching parameter and screening 
length are discussed in section 5 and 6, respectively. In section 7, we study the free en- 
ergy of the background black hole solution and discuss the phase transition between small 
black hole and large black hole in small chemical potential region. Section 8 is devoted to 
conclusions and discussions. 


2 Einstein-Maxwell-Dilaton system 


In this section, we use the potential reconstruction approach [12, 15] to study a 5D Einstein- 
Maxwell-Dilaton (EMD) system. In [15], the authors did not consider the coupling between 


gauge field and dilaton field in Einstein frame. Here we take the coupling into consideration 
in a more generic version 
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where the action (2.1) is written in string frame, Fuy = ô „Ay — 0, A, is the Maxwell field, 
Z(@) is an arbitrary function of dilaton field ¢ and Vg(@) is the dilaton potential. In 
Einstein frame we can rewrite the action as [12] 
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The Einstein equations from the action (2.2) read 
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where Ev = Ryv— sRow is Einstein tensor. We here consider the ansatz A = Ao(z)dt, ¢= 
o(z) for matter fields and 
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for the metric in string frame, where 7 = 1,2,3, Z is the radius of AdS; space, and A, is 


zZ 


the warped factor, a function of coordinate z. The metric in the string frame will be used 
to calculate the loop operator below. In Einstein frame the metric reads 
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with A, = A, — 26/3. In the metric (2.6), the (t,t), (z,z) and (zi, £i) components of 
Einstein equations are respectively 
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where b(z) = @?e?4¢/z?, and Ao(z) is electrical potential of Maxwell field. From those 
three equations one can obtain following two equations which do not contain the dilaton 
potential Vg(¢), 
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Eq.(2.8) is our starting point to find exact solutions of the system. Note that Eq.(2.8) 
in the EMD system is the same as the one in the Ejinstein-dilaton system considered in 
[13][14] and the last term in Eq.(2.9) is an additional contribution from electrical field. In 
addition, the equation of motion (EOM) of the dilaton field is given by 
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From equations of motion, once A,(z) is given, we can obtain a general solution to the 
system, which takes the following form 
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where ġo, Ago, Aoi; fo, fı are all integration constants and can be determined by suitable 
UV and IR boundary conditions. When Z(¢) = 1, the general solution reduces to the one 
given in [15]. Thus we have given a generic formulism to generate a set of exact solutions 
of the EMD system with a given A,(z). 


3 General asymptotical AdS black hole solutions 


Since here we are only interested in the black hole solutions with asymptotic AdS boundary, 
we impose the boundary condition f(0) = 1 at the AdS boundary z = 0, and require 
0(z), f(z), Ao(z) to be regular at black hole horizon zp and AdS boundary z = 0. There 
is an additional condition Ag(zp,) = 0, which corresponds to the physical requirement that 
Aga’ = g” Ay Ap must be finite at z = zp. 

We can parameterize the function f(z) in Eq.(2.12) as 
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We expand the gauge field near the AdS boundary to relate the two integration constants 
to chemical potential and charge density, respectively, 
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The temperature of the black hole can be determined through the function f(z) in 
(3.1) as 
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Following the standard Bekenstein-Hawking entropy formula, from the geometry given in 
Eq. (2.6), we obtain the black hole entropy density S as 
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where V3 is the volume of the black hole spatial directions spanned by coordinates x; in 
(2.6). 
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3.1 An analytical black hole solution 


In this subsection, we list an analytical solution of the Einstein-Maxwell-Dilaton system 
by using Eq.(2.12-2.15) with Z(¢) = 1. We impose the constrain f(0) = 1, and require 
o(z), f(z) to be regular at z = 0, and zp. We give the solution in Einstein frame as 
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where zo is an integration constant and V11, Vig are two constants from the dilaton potential 
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The two integration constants Vj; and Vi2 then can be expressed in terms of horizon zp 
and chemical potential u as 
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We can obtain the temperature of the black hole by using Eq.(3.6) based on the above 
formulas. In Figure [1], we show the temperature as a function of horizon radius zp, in 
cases of three different chemical potentials u. In this plot we take parameters ¢ = 1, zp = 
1,9, = 1. We see from Figure [1] that the temperature with respect to horizon zp is 
monotonic for a fixed chemical potential. A vanishing temperature means that the black 
hole is extremal with a smallest horizon radius. The smallest horizon radius increases as 
chemical potential becomes large. 


3.2 The hQCD model 


Based on the general solutions, in Ref. [15], a holographic QCD model is proposed to realize 
the confinement /deconfinement phase transition of QCD. Since the aim of this paper is to 


Figure 1. The temperature as a function of horizon radius z, for the analytical black hole solution 
with parameters £ = 1, zo = 1,9, = 1. 


further study the hQCD model, we here briefly review the main features of the model. In 
our hQCD model, the warped factor A;(z) takes the form 


A(z) = kz", (3.12) 
where k is a constant. To take this factor has various phenomenological motivations in 
order to build a successful holographic QCD model, for details see [15]. In this paper, we 
set k = 0.3GeV, following Ref. [15]. With this factor, we have the dilaton field ¢ as 

32,2 
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where we have set the integration constant ¢9 = 0, and H(z) is given by 
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The characteristic function of the black hole background takes the form 


ie tlie on L 


493" \ fo" gly) 3 dy o” g(a)dx 
0 g(x)dz 
_ 10 REI 3.15 
o” g(x)dz oe 
where 
gla) = reU) -3a (3.16) 


One can clearly see that the second term in (3.15) comes from the contribution of electric 
field. If one turns off the electric field, one can reproduce the black hole solution in Einstein- 
dilaton system [13]. In addition, the electric field A;(z) is given by 
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Figure 2. Plot (a): The black hole temperature as a function of horizon zp with different chemical 
potentials. When u > ue the temperature monotonically decreases to zero with increase of zp; when 
0 < u < He, the temperature decreases to a minimum at zm and grows up to a maximum at z),y and 
then decreases to zero monotonically. When p = He, one has zm = zm. Plot (b): The temperature 


of the black hole with u = 0.1GeV. The three black hole solutions with horizon z,, zp and z; have 


the same temperature. The black hole with zm < zp < zm is thermodynamically unstable. Here 
we take ggl = 1,k = 0.3GeV. In this hQCD model, we always fix k = 0.3GeV and accordingly the 
critical chemical potential is pe = 0.34GeV, which corresponds to the case zm = zy. 


The temperature behavior of the black hole solution is discussed in some details in [15], 
with respect to horizon radius and chemical potential. To be complete and for later use, we 
here briefly repeat some main salient features. In Figure [2] we plot the temperature with 
respect to horizon radius z, with different chemical potentials. One can see clearly that 
the temperature behavior crucially depends on the value of chemical potential: there is a 
critical chemical potential ue, beyond which the black hole is always thermodynamically 
stable, while when the chemical potential is less than the critical one, there is a region 
of horizon radius, where the black hole is thermodynamically unstable with negative heat 
capacity. To be more clear we plot in Figure [2](b) the temperature versus the horizon zp in 
the case u = 0.1GeV < uc as an example. One can see from the figure that the black hole is 
thermodynamically unstable in the region zm < zp < zm (branch 2 in plot b), where zm and 
zm are the black hole horizons corresponding to the minimal and maximal temperatures 
(Tm and Tm), respectively. In this region, the heat capacity of the black hole is negative. 
The black hole solutions in the regions zp < Zm and zp, > zm are thermodynamically stable 
(branch 1 and branch 3). When p > He, 2m and zm are degenerated to one point. Note 
that in contrast to the case in Figure [1], there are local minimal and maximal values of 
temperature in small u cases. This is crucial to realize the critical point in the T — u phase 
diagram of the hQCD model [15]. 

At this stage let us mention that it is quite interesting to compare the temperature 
behavior of our black hole solution (see Figure [2]) with the one for the Reissner-Nordstrém- 
AdS black hole (see Figure 3 of Ref. [52]). They look quite similar. But there exist 


essentially different points between them. At first, in our case, the temperature is plotted 
with respect to horizon radius zp, for a fixed chemical potential, while in [52] it is plotted 
with respect to horizon radius for a fixed charge (note that there the horizon radius r+ 
corresponds to the inverse of horizon radius zp, here). Second, our black hole solution has 
a Ricci flat horizon, while the black hole horizon discussed in [52] is a sphere. The third is 
that we will discuss the phase transition in grand canonical ensemble (for a fixed chemical 
potential) while the phase transition between a small black hole and a large black hole 
discussed in [52] is in canonical ensemble (for a fixed charge). As we will see shortly that 
just due to the similarity of temperature behavior, in our case, the phase transition will 
also appear between a small black hole and a large black hole in grand canonical ensemble 
when the chemical potential is less than the critical value. 


4 Drag force 


In this section, mainly following [18] we compute the drag force experienced by an ex- 
ternal probe quark traversing in the QGP in the present hQCD model. The holographic 
computation of the drag force is based on the gauge/gravity correspondence between the 
deconfined phase of hQCD model at finite temperature and its dual garvity realized as 
an aAdS black hole solution in the Einstein-Maxwell-Dilaton background. In this aAdS 
black hole geometry, the gravity dual of the probe quark is described by an infinitely long 
fundamental string. One of its ends is attached to the boundary of the bulk spacetime. 
The body of the string extends along the radial direction and the free end of the string 
goes parallel to black hole horizon. The gauge/gravity duality suggests an identification 
between the end point of the string and the probe quark. Furthermore, the body of the 
string captures the effects of thermal plasma through which the external quark is mov- 
ing. In this dual gravity picture, the string trails back and imparts a drag force on it’s 
endpoint that is attached to the boundary. This drag force is obtained by calculating 
the rate of change in string momentum. The drag force is a function of temperature and 
chemical potential. The boundary gauge theory we are considering is on M4 described by 


2 z. The dynamics of a fundamental string is completely 


the boundary coordinates t, x1, x 
specified by the Nambu-Goto action in the black hole background within the string frame 
(2.5). On this background the world sheet action reads 
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where gag is the induced metric on the world sheet and Gy is the background metric. The 
equation of motion derived from (4.1) is given by 
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where Aa is the covariant derivative with respect to gag and P% is the world sheet current 
of space time energy-momentum of the test string. We consider the motion of the string 


along x!. In the gauge, T = t and o = z, the string dynamics can be completely specified 
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by the function x!(t, z). In this case, the Lagrangian reads 
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To capture the dragged motion of the quark in the boundary theory we assume the following 
ansatz in the bulk [18] 
a(t, z) = ut + €(z). (4.5) 


Here we have assumed only the late time behavior of the string motion. With this ansatz 
the Lagrangian reduces to 
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The momentum which conjugates to €(z) reads 
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For the sake of consistency it is important to invert the equation (4.7) and write it in the 


following way 


(4.8) 


Here the positive sign is taken due to the trailing nature of the string profile [18]. To obtain 
the string profile we have to solve the differential equation (4.8). To have a real €(z), we 
further impose the constraints 


F) lz=z = f, 
Me? |p», = 


| 
e 


EPEL (4.9) 
so that one has €’|,-,, = v?/f?, keeping finite. The profile of the string is defined in the 
region with z < z,, That is, there is a maximal value z, < zp for the string profile. 

The constraints are very useful to figure out the final form of the drag force. Before 
to compute the drag force, we here mention the relation between the drag force in the 
boundary field theory and the dissipation of momentum flowing down the string, in light 
of AdS/CFT correspondence. In the boundary theory the presence of the thermal medium 
results into dissipation of energy and momentum of external quark until it reaches thermal 
equilibrium with the medium. In the bulk theory the momentum is flowing down the string 
from the boundary to the bulk and the change of momentum at a given spatial point on the 
world sheet for a given time interval can be calculated. The identifications of the endpoint 


=] i= 


of the string attached to the boundary with the quark and of the string in the bulk with 
the thermal medium around the quark suggest that the drag force can be realized in terms 
of the force imparted by the string on its boundary endpoint. To calculate the change of 
string momentum due to its motion along xı direction, we consider a closed curve on the 
world sheet and study how the momentum is conserved around this curve [53]. According 
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Figure 3. This plot shows a closed path in an anti-clockwise direction on a world sheet bounded 
by coordinates [A = 7), 02], [B = 72,02], [D = 72,01] and [E = 7,04]. 


to the conservation of world sheet current of space time energy-momentum of the test 
string, the total flux calculated around the path C must be zero, 


(Piao — Pī dr) = 0. (4.10) 


Note that one end of the string is attached to the boundary and the other end close to 
horizon is free. In the static gauge, Eq. (4.10) reduces to 
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where Pha is the zı component of the total momentum at time t. Consequently the drag 


force is defined as 
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Finally we have to replace all the gravity parameters in terms of gauge theory parameters. 
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Before doing that, we analyze the form of the constraints case by case. The exact forms of 
constraint for the solutions (3.9) and (3.15) are given respectively by 
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It is always desirable to express the drag force in closed analytic form as a function of 
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gauge theoretical variables. However it is very difficult to obtain analytic forms for the 
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Figure 4. This plot shows the drag force from the analytic black hole solution as a function of T 
for chemical potential u = 0.01,0.80, and 1.20 respectively. Here we take v = 0.1. 
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Figure 5. This figure shows the drag force as a function of T for the chemical potential u = 0.10 
(plot a), 0.34 (plot b), and 0.80 (plot c), respectively, in the hQCD model. Here the dashed curves 
stand for the behavior of drag force in confined phase which is denoted by CF in the figure. In the 
confined phase, in fact the drag force is not well defined, meaning that the dashed curves do not 
make any sense here. 


constraint (4.13) and the temperature (3.6). Instead we here solve them numerically and 
plot the drag force with respect to gauge theory parameters, e.g, temperature and chemical 
potential, so that the qualitative features of the drag force can be revealed. 

Certainly the analytic black hole solution (3.9) does not dual to a QCD model. As 
a warm-up exercise, we plot the drag force in Figure [4] for a model dual to this black 
hole solution. We see that the drag force monotonically decreases with temperature and 
for a fixed temperature it becomes large as chemical potential grows. These features are 
qualitatively expected in realistic QCD. However, the behaviors of jet quenching parameter 
and screening length in the solution (3.9) are far away from QCD phenomenon and therefore 
we do not consider this solution from now on. 

In Figure [5] we plot the drag force for our hQCD model given by the solution (3.15) 
with different chemical potentials u = 0.10, 0.34, and 0.80, respectively. We can see from 
the figure that for fixed chemical potential and temperature, the drag force increases with 
the velocity of the quark, while for fixed chemical potential and velocity, the drag force 
increases with temperature. These are expected features in QCD. In particular, let us 
note that in the low temperature region with small chemical potential, the drag force is a 
multi-valued function of temperature [see plot (a) and (b)|, while it becomes a monotonic 
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Figure 6. The figure shows the drag force as a function of T for three chemical potential u = 
0.10, 0.34, and 0.80, respectively in the hQCD model. The dashed parts of the curves stand for the 
drag force in the confined phase which denoted by CF in the figure. Here we take v = 0.1. 


function with large chemical potential [see plot (c)]. This feature is closely related to the 
confinement /deconfinement phase transition in this hQCD model [15]. The dashed parts 
of curves in plot (a) and (b) denote the drag force in the confined phase and actually they 
do not make any sense here since drag force is not well-defined in the confined phase. Our 
result for the drag force in the deconfined phase is in agreement with the one in [54]. For 
comparison, in Figure [6] we plot the drag force versus temperature with three different 
chemical potentials u = 0.10, 0.34 and 0.80, respectively. In this figure the velocity of 
quark is taken as v = 0.1. 


5 Jet Quenching parameter 


In this section, we use the AdS/CFT duality to compute the jet quenching parameter in 
our hQCD model. The holographic method to compute this quantity needs a consideration 
of Wilson loop (C) traced out by a qq pair [33]. The Wilson loop is taken to lie along the 
light cone in the gauge theory. The gravity dual of this gq pair is represented as the two end 
points of a fundamental string, attached to the boundary of the bulk spacetime. Body of 
the string hangs along the radial direction and up to the horizon of aAdS black hole. The 
Wilson loop is mapped, in the dual theory, as the string world sheet. The jet quenching 
parameter (ĝ) is related to the thermal expectation value of the light-like Wilson loop op- 
erator, (W(Cyight—like)) [55]. The holographic correspondence between thermal expectation 
value of the light-like Wilson loop operator in fundamental representation (WF (Clight-like)) 
and the exponential of the worldsheet action, e75, leads us to obtain a working formula 
of ĝ in dual gravity theory. Here, S stands for the worldsheet action of the fundamental 
string. 


(WF (Ctight—tike)) = exp[—S(C)]. (5.1) 
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In the planar limit, considering the fact Tr aq; = Tr? na » the relation between W‘ (C) and 
W4(C) can be easily established as 


(WA(C)) = (WFC))’. (5.2) 


Now we start with the background black hole solution in string frame 


2 
A + da'da) + dx*dx? + dažda”) : (5.3) 
z 


By introducing the light cone coordinates defined as 


— (roa + 


the black hole metric (5.3) can be rewritten as 


D2 As Pi z 
d3 = £ (Eart (de)”) — (1 + f(2))(dr* dx) 


z2 2 
dz? 27.2 34.3 
+— + da*dx* + dx’dx”). (5.5) 
f(z) 
We take the gauge with T =27(0< £7 < L7), o =a?(-@ < z < 4B), and set the pair 
of quarks at z? = +2 on zt = constant, x? = constant plane. In the limit with L~ > Lə 


the string profile is completely specified by z = z(c). Following [33], using (5.1), (5.2), one 
has 


(wA(c)) = exp(—_ gal £3), (5.6) 


where the jet quenching parameter is defined as 


1=— e ° 
where S' is the Nambu-Goto action of the string and So is the self energy from the mass of 
two quarks. 

Substituting the induced metric of the fundamental string into the Nambu-Goto action 


(4.1), we get 


1 
S= -z7 | drdo/=aetgan, 


£2 12 
Le f 2 e2As z 
= do 1— f(z) + =). 5.8 
Fao f° vE Eo Fp) (5.8) 
Since the integrand in (5.8) does not explicitly depend on ø, one can regard o as time and 
the integrand as a Lagrangian. In this case the corresponding Hamiltonian is conserved. 
That is, we can have 


ƏL, 
ay? LE, (5.9) 
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where F is a constant and £ is the integrand in (5.8). From this relation we obtain the 
equation of motion for z as 


ee (= f(z) 


2# = FA ie. (5.10) 


We choose the boundary conditions z(+2) = 0 and z'(0) = 0. In that case, the 
turning point zp is determined by solving Eq.(5.10). Since z'(c) is a real function, so the 
square of it should be non-negative. The realization of boundary condition z (0) = 0 at 
the turning point requires the proper choices of zeros and the positivity region of the right 


hand side of Eq.(5.10). From the boundary conditions of the black hole solution 


hin (2) =0,. lim f(z) =1, (5.11) 
22, z—0 
together with the fact that we are interested in the case with small Ẹ, it is clear that the 
factor “0-sp-#"-" is always positive near the black hole horizon and negative near the 


2 


boundary. To remove the region with a negative z2, we consider a modified boundary at 


z= ô. We assume that at z = Zmin, 


e4As(2min) (1 — Flamin) ) 


min 


—~1=0, (5.12) 


4s If) _ yy 
a-a) 


and ô > Zmin. In the region 6 < z < zp, thus, the factor [=> is always positive. 


So only viable solution of z? = 0 is 
fe) =0> r= 2). (5.13) 


That is, the turning point is just at the horizon. The distance between two quarks can be 


L f” E 
= 1 Tama fet By Ft (5.14) 


As we are interested in the small Lə limit, considering the smallness of Æ, we can expand 
Eq.(5.14) in terms of E as 


determined by 


L Zh 2 ,—2As E? Zh —6A, _6 
= go | ga On: (5.15) 
w Is “Fa 2 Is “FOP 

Inverting (5.15) suitably, we can obtain E up to the leading order of Lə as 


Lo 


B= e a (5.16 
4 f(-f) 
Thus we can obtain the string action 
Loe fr et^ (1 — f) 


S (5.17) 
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Clearly this action is divergent. The divergence comes from the contribution of mass of 


two quarks. With the gauge x~ = 7 and z = øg, the self energy of two free quarks reads 
L- 2 Zh 2As T= 
2 = | dz w= (5.18) 
Vra 5 Z Í 


Thus the regularized action up to the leading order of Lə is given by 


L- LE 1 4 
2 a FO (5.19) 
8 ora! Zh d z2e—2As 
vore fz" dz Tf 


Sr = S — So = 


With the definition of the jet quenching parameter (5.7), we finally reach 


0 e 1 
j f(1-f) 
In fact the cutoff here can be removed by noting the fact that the integrand is regular 


inside the region 0 < z < zp, i.e, from the horizon to the real boundary, 


Zh 2 5—2As Zh 2 5—2As ô 2 4—2As 
AE ze Ze 
f d (5.21) 
ô 


z—— = dz ———— — dz ———. 
Vif) Jo vfA-f) Jo vVIQ-P) 


The second integral in the right hand side of the above equation smoothly vanishes in the 
limit 6 > 0. So the final expression for the jet-quenching parameter is 


g 1 
— 7 2e-2A3 ` (5.22) 


Because the black hole metric is still too complicated to obtain an analytical expression of 
the jet-quenching parameter in terms of physical parameters, we plot in Figure [7] the jet- 
quenching parameter as a function of temperature in the hQCD model with three chemical 
potentials u = 0.10, 0.34 and 0.80, respectively. For large u > ue cases, the jet-quenching 
parameter decreases monotonically with temperature, which agrees with the one in [56] 
qualitatively. On the other hand, when u < ue, the jet-quenching parameter is a multi- 
valued function of temperature in low temperature region and it decreases monotonically 
with respect to temperature in high temperature region. The multi-valued behavior of the 
jet-quenching parameter in low temperature region is clearly related to the first order phase 
transition between hadron phase (confined phase) and QGP phase (deconfined phase). The 
jet-quenching parameter confirms the hydrodynamical description of QGP phase and agrees 
with the real QCD expectation in high temperature. Once again, as the drag force in the 
confined phase, the dashed parts of curves in Figure [7] denote the jet-quenching parameter 
in the confined phase and thus they do not make any sense. 


6 Hot plasma wind and screening length 


The screening length is defined as the maximum length achieved by a quark-antiquark 
bound state at temperature T > Te, beyond which the pair dissociates. For quark-antiquark 
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Figure 7. The figure shows the jet-quenching parameter as a function of T for three chemical 
potentials u = 0.10,0.34, and 0.8, respectively, in the hQCD model. The dashed parts of curves 
stand for the jet-quenching parameter in the confined phase which are denoted by CF in the figure. 


pair, the energetically favorable configuration in the dual gravity theory is a fundamental 
string with both ends attached to the boundary. The attached endpoints correspond to 
the qq pair whereas being separated beyond the screening length, thus dissociated from 
each other, the pair maps into two separate strings hanging from the boundary. In [53], 
the screening length is computed in the rest frame of gq pair and the plasma wind flows 
at a constant speed v for the hot M = 4 SYM plasma. This setup is identified with a 
quark-antiquark pair moving in hot M = 4 SYM plasma. In this section, we compute the 
screening length for the hQCD model (3.15) in the same way as in [53]. 

In the static frame of gq pair, we assume that the hot plasma is moving with velocity 
v in the negative «3 direction. The Wilson loop we are interested in lies in the t— xı plane 
specified by the length 7 and L respectively. We assume 7 >> L such that the string world 
sheet is invariant under translation along the time direction. The boost we are considering 
is defined as 


dt 
dzz = — sinh ndt + cosh ndrs, (6.1) 


cosh ndt — sinh nds, 


where coshn = y, sinhy = yu and y = 1/1 — v? is the Lorentz boost factor. With the 
Lorentz transformation, we obtain the boosted black hole metric in string frame 
ds}, = H(z)[-(1 — (1 — f) cosh? n)dt? + (1 + (1 — f) sinh? n) (dz), 
2. dz? 


—2(1 — f) cosh n sinh ndtdx? + (dx)? + (dz?) + OL (6.2) 
where H(z) = ¢?e?4s /z?. We prefer to work in the static gauge 
T = t,o =2',2*(c) = x°(c) = constant, (6.3) 
with the following boundary conditions 
L , 
z(o +5) 0, z(o = 0) = zaz (o =0) = 0. (6.4) 
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Thus the world sheet metric induced on the boosted background is given as 


Orr = —H(z)(1 — (1 — f) cosh? n), 


Jro = Jor = 0, 
a 
Goo = A(z) 1 + (1+7 6.5 
(z)[1 + ( T z (6.5) 
Then the Nambu-Goto action for the string takes the form as 

T 5 12 
S=-— | doH(z)\/(1- (1 f) cosh? n)(1 + =). (6.6) 

TQ 0 f 


As the Lagrangian £ in (6.6) does not depend on o explicitly, the corresponding Hamilto- 
nian is conserved and can be viewed as a constant of motion 


OL , 


With this we can cast the equation of motion in the form as 


fP- (1 = f)cosh”n) =g] 
2 = 1 (6.8) 
q 

It is evident from the constraint (6.8) that at the horizon, z = zn, where f(z,) = 0, the 
factor al = (1 — f) cosh? n) 1 = E sinh? 7 — 1 is always negative. At the boundary, 
f (0) = 1, the factor Ba — (1 — f) cosh? n) — 1 = +- — 1 is always positive for small values 
of q < H. Therefore in the range 0 < z < zp there must be a location (z = ze) where 
Ba — (1 — f) cosh? n) — 1 switches its sign. Accordingly z = ze is the physical turning 
point of the string configuration. The string can not be stretched up to the horizon as 2’ 
is an imaginary quantity in the region ze < z < zp. By solving the equation 


f (Zc) H? (ze) cosh? n E H? (ze) sinh? n 7 


1=0, 6.9 
2 7 (6.9) 


the turning point can be numerically determined. Then one can obtain the binding energy 
between the quark and antiquark pair through calculating the action (6.6) with constraint 
(6.8) 


(6.10) 


where So is given by 


Sod deal T (6.11) 


Ta Jo 
The distance between quark and antiquark can be calculated from (6.8) as 
1 


L Á 
E a S (6.12) 
2q / H4/ f[(1 — cosh? n(1 — f)) — red 
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Figure 8. Plot (a) shows the quark-antiquark distance as a function of q for a fixed rapidity, while 
plot (b) shows the binding energy with respect to the distance. In both plots we fix the chemical 
potential u = 0.10, 0.34, and 0.80, respectively. We have set a same temperature T to obtain these 
curves. 
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Figure 9. The screening length versus the rapidity 7 for the cases u = 0.1, 0.34, and 0.8, respec- 
tively. We have set a same temperature T to obtain these curves. 


It is not possible to work out the integration in (6.12) explicitly. To determine the screening 
length, we plot the distance L with respect to the constant of motion q for a fixed rapidity 
n in Figure [8] (see plot (a)). It turns out that for a fixed value of rapidity, there exists 
a maximum for L, which is regarded as the screening length Ls = Lmazx(n)/(aT). Plot 
(b) in Figure [8] shows the binding energy V given by (6.10) with respect to L. One can 
see from plot (a) that the quark-antiquark distance starts from zero when q is also zero, it 
increases sharply with respect to q, reaches its maximum at a certain q, and then decreases 
monotonically to zero at some finite q. In between these two zeros, there exists a single 
L = Lmax beyond which there is no solution of Eq.(6.12). This implies the quark-antiquark 
pair dissociates beyond L = Lmar. We identify Lmaz(n)/(nT) with the screening length 
Ls. For the u = 0.1 case, Lmaz ~ 1.4 and Ls ~ 1.4/(nT) œ 0.45/T, close to the lattice 
calculation Ls ~ 0.5/T [57] of the static potential between heavy quark and antiquark in 
QCD. Plot (b) shows that there are two branches for the binding energy in the region 
L < Imax. The branch with dashed curves has a higher energy than the one with solid 
curves. This implies that the branch with dashed curves is physically disfavored. 

The screening length L,(7) as a function of rapidity can be obtained numerically as 
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illustrated in Figure [9]. One finds that it decreases with velocity which indicates that the 
quark-antiquark pair dissociates at a lower temperature as it is moving. This behavior is 
also observed in [54]. If the qualitative behavior holds for QCD, it will have the consequence 
for quarkonium suppression in heavy ion collision. Additionally, our results show that the 
case with smaller chemical potential has a larger screening length. 


7 Free energy and phase transition 


In this section, we would like to calculate the free energy of the black hole solution dual to 
the hQCD model by following [46]. Ref. [46] studies the free energy of a generic Einstein- 
Maxwell-Dilaton system. The system considered in [46] is given by 


T= Mp" fy dt a/—g |R- EPF? — 300) + V(®)| + Ton, 


(7.1) 
logg = 2M9! Jan d’zy -hK , 
with following ansatz 
2 2A(u) z 2 igni du? 
ds? =e (—F(wat + dride ) BT A=A,(u)dt, @=O(u), (7.2) 
u 


where Mp is the Planck mass, and K is the extrinsic curvature of the finite boundary 0M 
with induced metric h. By the following transformations 


d -iui # = edwin (L ay Flu(z)) = fle), 


du dz’? du? dz? dz dz 
auj = Syla), Z(6)=1, d=4, Mi= wae 
Haas Vea, (7.3) 


we can change the system and ansatz to ours discussed in this paper. Here we have chosen 
A(u(z)) = Ae(z) — log(z) as a gauge. Varying the action (7.1) yields equations of motion 
for the gravitational field 

1 Z(®) 


1 1 1 


The Maxwell equation and the explicit forms of Einstein equations read 


d zo. 
ia (d—2)A — 
z- (e ZA,) 0, (7.5a) 
2(d-1)Ã +° =0, (7.5b) 
f+dAf —e?4Z A? = 0, (7.50) 
DS 2 l a 1 es 
(d—1)Af + (aa A, = 5%) f-V+ gee At = IN (7.5d) 


where an over dot stands for the derivative with respect to u. Eq.(7.5a) is the Maxwell 
equation and the other three Eq.(7.5b) (7.5c) (7.5d) are obtained from Einstein equations. 


es e 


One can easily check that Eq.(7.5b) and Eq.(7.5c) correspond to Eq.(2.8) and Eq.(2.9), 
respectively. In addition, the equation of motion for the scalar field is 


1 / Z2(@) 29 
7z (/—g0u®) + V'(®) — =o =0, (7.6) 


where the prime denotes the derivative with respect to ®. From (7.5a), one obtains 


u 


p : 
where u and p are the chemical potential and charge density of the black hole solution, 
respectively, and ug stands for the UV boundary. 

By defining the superpotential W in the following way [46] 


= W'(®), (7.8) 
the equation (7.5b) can be solved as 
A W(®) 
A =-=. . 
2(d — 1) (a 
Equivalently, 
: > 1 ? . W(e) ? 66 
A(®) = i-z). ete | 7.10 
(®) = Ao -T Ja È WO) o+ a WO (7.10) 


where Ag = A(&p), (u = ug) = Bo. 
The temperature T associated with the black hole in metric (7.2) is given by 


(7.11) 


U=Up, 


= 1 re 
T=—|e4 
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where up denotes the black hole horizon. The entropy density for the black hole solution 
(7.2) is 


Š = Aneli -DÂ@r) — gre 40-4 Jas W, (B, := B(up)). (7.12) 


where Met = ] has been taken. 

The on-shell action J has a bulk contribution Ig and a boundary contribution IGH 
from the Gibbons-Hawking term. Here we evaluate it with an UV cutoff at u = uo or 
z =e. A standard computation gives us with the action density 


Joe shell = T1(-W T PyetAuo)) u. (7.13) 


With the equations (7.3)(7.5c) and A;(up,) = 0, the free energy density is found to be 
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Figure 10. The free energy of the background black hole solution with respect to temperature. It 
shows that in the small u < ue case, the free energy has a multi-valued behavior, while it is absent 
in the large u > pic case. The parameters l = 1, 167G's5 = 1, gg = 1, k = 0.3GeV are taken here. 


where b(z) = e4e()/z. Note that in the second line of Eq.(7.14), we have considered the 
case d = 4 and the expression of the free energy is written in the ansatz (2.6) in Einstein 
frame. The first term in Eq.(7.14) corresponds to the black hole mass and the last term 
can be expressed as — up. 

With the analytical result (7.14), one can obtain the free energy of the black hole 
solution (3.15) with respect to temperature with various chemical potentials numerically. 
Figure [10] shows the free energy density of the black hole solution. It can be seen that 
the free energy density is a multi-valued function of temperature in the case u < ue and it 
becomes monotonically in the large u region u > ue. The multi-value behavior of the free 
energy indicates the existence of a first order phase transition in the small u region. 

To clearly see this, it should be helpful to recall the temperature behavior of the 
black hole (see Figure [2]). One can see that when u > fc, the black hole is always 
stable with positive heat capacity. When u < He, one can see from plot (b) that as 
T > Tw, the black hole is also always stable with positive heat capacity (in branch 1), 
when 0 < T < Tm, the black hole is also stable with positive heat capacity (in branch 3), 
while when Tm < T < Tm, there exist three black hole solutions with a same temperature. 
For instance, for temperature Tp, there exist three black hole solutions with horizon radii 
Zp, % and zp, in branch 1, 2 and 3, respectively. The black hole solutions in branch 2 
are unstable because those black holes have negative heat capacity. We call black holes in 
branch 1 and branch 3 as large and small black holes, respectively, when Tm < T < Ty, 
because the black holes in branch 1 have large horizon radius and entropy, while the black 
holes in branch 3 have small horizon radius and entropy. From Figure [2,10], we can see 
that when u > ue, the black hole solution always dominates, which indicates the dual field 
theory is in deconfined phase, while in the u < pe regime, when Tm < T < Tym, there 
exists a phase transition between large black hole and small black hole. In this transition, 
the entropy has a jump, which means this is a first order phase transition. On the other 
hand, when 0 < T < Tm, the small black hole solution in branch 3 is dominated. This 
is quite different from the case of Schwarzschild black hole in AdS space. There when 
temperature is less than a certain value, there does not exist black hole solution. In this 
case, the thermal gas is dominated. Above a certain temperature, there exist two black hole 
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solutions, one is stable and the other is unstable. Therefore there exists a well-known first 
order phase transition (Hawking-Page phase transition) between black hole and thermal 
gas in AdS space [47]. This transition is explained as the confinement /deconfinement 
phase transition of dual field theory. In our case even in the low temperature regime 
with T < Tm, the small black hole is dominated, the free energy for these small black 
holes is negative. Compared to thermal gas, those black hole solutions are still favored. 
When temperature increases beyond Tm, the black hole solution with large horizon radius 
in branch 1 is favored than the one with small horizon radius in branch 3. In this case, 
a phase transition happens. Since the entropy of small black holes in branch 3 is much 
small than the one of large black holes in branch 1, this phase transition therefore can be 
identified as the confinement /deconfinement phase transition in the dual theory. Of course, 
this phase transition is not same as the Hawking-Page phase transition. The free energy 
calculation further confirms this picture. And this picture is completely consistent with 
the one in [15] by computing heavy quark potential in this model. In addition, we would 
like to stress that the phase transition between large black hole and small black hole is the 
same as the one found in [52], but the transition in our case is in grand canonical ensemble, 
while the one in [52] is in canonical ensemble. As a result, T — Q phase diagram in Figure 
1 of [52] is similar to our T — u phase diagram. 

In a word, in the small u region there exists a first order phase transition, while it is 
absent in the large u region through studying the free energy of the system. The existence 
of the first order phase transition in the small u region is consistent with the analysis 
of Wilson loop in this black hole background [15]. The existence of the critical point at 
U = Uc is in agreement with recent lattice calculation given in [58]. Let us notice that in 
[15] we claim that there is a crossover or higher order phase transition in the large u > pe 
region, while the free energy calculation here shows no such phase transition. In this sense, 
free energy may show less information about phase transition and one should take more 
probes into considerations, such as Polyakov loop [15], to understand the phase structure 
completely. These two conclusions are in fact not in contradiction with each other. The 
reason is as follows. The free energy calculation here only concerns with the black hole 
background, from the point of view of QCD, the black hole background does not include 
the degrees of freedom of quarks, while the Wilson loop calculation in [15] is related to the 
dynamics of quarks. 

In Figure [11] we plot the phase diagram of the hQCD model. The difference between 
the confinement /deconfiement phase transition lines come from the different consideration 
as mentioned above. 

At the end of this section, let us stress the question whether there exists the so-called 
Hawking-Page phase transition between black hole and thermal gas in our model. Note 
that in the above analysis, our phase transition happens between small black hole and large 
black hole in the region with small chemical potential u < ue, such a transition does not 
occur as > He. It is therefore of some interest to investigate whether there exists any 
transition between black hole solution and thermal gas solution in our hQCD model. For 
this aim, we have to first find the thermal gas solution. Unfortunately, in this potential 
reconstruct approach, to find out the thermal gas solution is not an easy job since we 
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T/GeV 


Figure 11. The phase diagram of the hQCD model. The red and dashed red curves are obtained 
in [15] through studying VEV of Polyakov loop. The two curves stand for first order and continuous 
phase (or second order) transition respectively. The thick black and blue curves correspond to the 
minimal temperature Tm and maximal temperature Tm in small u region. The thick black curve 
is obtained by calculating the free energy. The black dot denotes the critical point. T, denotes the 
transition temperature. 


have not an explicitly analytic form of the potential for the dilaton field in our hQCD 
model. But even so, we still can make some arguments which lead to the conclusion that 
there does not exist any phase transition between black hole and thermal gas solutions 
in the hQCD model. The main reason for this comes from the fact that in the hQCD 
model, for any temperature, there exist corresponding black hole solutions. This can be 
seen from Figure [2]. It shows in plot (a) that when u > ue, the black hole temperature 
increases monotonically from zero with the smallest black hole radius, which corresponds 
to an extremal black hole, while when u < ue, plot (b) shows that as T < Tm, there exists 
a region with stable black holes with small horizon radius, including the extremal black 
hole; and as T > Tm, there exists the region with stable black holes with larger horizon 
radius. In addition, for a given temperature, it is generally believed that the black hole 
solution has much large entropy than the thermal gas, so that the black hole phase is always 
dominated. As a result, whatever the chemical potential is, the black hole solutions are 
always dominated over the thermal gas solutions. This is quite different from the case of 
Schwarzschild-AdS black hole. In the latter case, there exists a minimal temperature, below 
which there is no black hole solution. This leads to the conclusion that in low temperature 
the thermal gas solution is favored, while in high temperature the black hole solution is 
dominated, and thus there must exist the Hawking-Page phase transition between the black 
hole solution and thermal gas solution at a certain temperature. Therefore in our case, black 
hole solutions are always dominated and there does not exist any phase transition between 
black hole solution and thermal gas solution. This might be a common feature for black 
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hole solutions with Ricci flat horizon in AdS space. But as we analyzed above, there exists 
indeed the phase transition between small black hole and large black hole in our model with 
u < He. This phase transition can also be understood as the confinement /deconfinement 
phase transition in QCD since clearly the small black hole has a much less entropy than 
the large one. And this interpretation is also consistent with the calculation of heavy quark 
potential in [15]. 


8 Conclusion and discussion 


In this paper, we have continued to study the holographic QCD model proposed in [15], 
in an Einstein-Maxwell-Dilaton system. At first we have generalized case with a non- 
minimal coupling between Maxwell field and dilaton field, and given a generic formulism 
for generating a set of exact and asymptotic AdS black hole solutions in the EMD system. 
After briefly reviewing the main features of the hQCD model, we have studied some aspects 
of QGP phase of the hQCD model by calculating some quantities such as drag force, jet 
quenching parameter and screening length. The calculations show that the behaviors of 
those quantities are consistent with the expectation from realistic QCD. 

It is found that the drag force increases monotonically with temperature which is quite 
good consistent with real QCD phenomenon in the larger chemical potential region with 
HU > He. In the small chemical region with u < pe, the drag force also monotonically 
increases in the high temperature region, while in the low temperature, it shows a multi- 
valued behavior. Note that in the case u < ue, the solution is dual to the confined phase 
of QCD. In that case, the drag force is not well defined. Therefore the change from 
the multi-valued behavior to the monotonic behavior just manifests the existence of the 
first order phase transition. The jet quenching parameter has monotonically decreasing 
behavior versus temperature, which is also consistent with QCD experiments in p > He. 
For the u < He case, the jet quenching parameter agrees with real QCD expectation in high 
temperature and once again, it shows the multi-valued behavior in the low-temperature 
region. As in the case of drag force, the multi-valued behavior of jet quenching parameter in 
the low temperature region is consistent with the existence of first order phase transition in 
this hQCD model. For the screening length we have plotted the separation between quark 
and anti-quark with respect to the constant of motion q. It is clear from the plot that for 
both cases when u > pie and u < He the dipole dissociates beyond a maximum separation 
distance, namely the screening length Ls. We have also calculated the binding energy as 
a function of separation distance. In addition, we have presented L;(7) and found that 
there are qualitative consequences for quarkonium suppression in heavy ion collisions in 
this hQCD model. 

We have calculated the free energy of the background black hole solution and further 
confirmed that there exists a first order phase transition in small u region between small 
black hole and large black hole [15]. When u > ue, the phase transition is absent from 
the point of view of the free energy. The existence of the critical point is consistent with 
the result in [15]. Further we have argued that there does not exist any phase transition 
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between black hole solutions and thermal gas solutions in this model, and the main reason 
is given. 

In this work we have studied some aspects of QGP phase in the hQCD model proposed 
n [15]. The results are encouraged and are consistent with the expectation of real QCD. 
Thus it would be of great interest to further investigate the hQCD model. For example, 
it is required in the model to study the spectra of hadrons, chiral phase transition and 
its phase diagram [59][60][61], hydrodynamical properties of QGP, and color flavor locked 
phase [62][63][64], etc. Furthermore, quantum corrections [65] to above physical quantities 
are also deserved to consider in the coming works. 

Finally we would like to mention that in this paper the drag force, jet quenching and 
screening length have been calculated through the Nambu-Goto action of fundamental 
string in the AdS/CFT correspondence. As one knows that the Nambu-Goto action makes 
sense only in string theory. Therefore it becomes a crucial issue whether our hQCD model 
can be embedded into some string theory. This is also a key point for a large kind of 
phenomenal models of holographic QCD. At the moment, we cannot show that the model 
could be a consistent truncation of some low energy effective theory of string theory, but 
we hope this model is helpful to build a holographic model dual to the realistic QCD. 
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